Spontaneous symmetry breaking in thermalization and anti-thermalization 
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The phenomenon of spontaneous symmetry breaking is investigated in the dynamic thermalization of a de- 
generate quantum system. A three-level system interacting with a heat bath is carefully studied to this end. It is 
shown that the three-level system with degenerate ground states might have different behaviors depending on the 
details of the interaction with the heat bath when the temperature approaches zero. If we introduce an external 
field to break the degeneracy of the ground states and let it approach zero after letting the temperature approach 
zero, then two possibilities will arise: the steady state is a definite one of the degenerate states independent of 
the initial state, or the steady state is dependent on the initial state in a complicated way. The first possibility 
corresponds to a spontaneous symmetry breaking of the system and the second one implies that the heat bath 
could not totally erase the initial information in certain cases. 
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Introduction. — Conventionally, thermalization is under- 
stood as a dynamic process in which an open quantum system 
(OQS) approaches an equilibrium state with the same temper- 
ature T as that of its heat bath [1]. According to the third 
law of thermodynamics (the generic version), the entropy of 
a non-degenerate system will vanish at the absolute zero tem- 
perature. That is to say, at the absolute zero temperature, the 
system will reach a steady pure state. Is it still the case when 
the system has degenerate ground states? This turns out to 
be a subtle problem related to the process of thermalization. 
In quantum physics, thermalization seems to be a much more 
complex concept than in classical physics. In fact, most re- 
cently a new kind of thermalization, called canonical thermal- 
ization UIH, has been proposed. 

At zero temperature a system with degenerate ground states 
might have a finite entropy depending on the degree of de- 
generacy d: S — kg\nd [6]. Roughly speaking, this can be 
understood by studying the following two non-commutative 
limit processes: letting the perturbation introduced to break 
the degeneracy approach zero and letting the temperature of 
the heat bath approach zero. Actually, in the thermal equi- 
librium case (the system is thermalized to reach a thermal 
equilibrium state), if we first take the second limit and then 
take the first limit, the OQS will reach a definite pure state 
and thus have a vanishing entropy. This is a phenomenon of 
spontaneous symmetry breaking (SSB) Q @]. On the other 
hand, if we reverse the order of these two limit processes, the 
OQS will reach a maximally mixing state of all the degenerate 
ground states and thus possess a non-vanishing entropy. 

Unfortunately, the above discussion about the thermaliza- 
tion of a system with degenerate ground states proves to be 
overly simplified. In the study of the dynamic thermalization 
of a simple system, we find that the above mentioned SSB in 
the thermalization can only happen when bath induced tran- 
sition is not forbidden by some selection rule. If there is a 
selection rule to forbid the bath induced transition between 
the degenerate ground states, the steady thermalized state will 
depend on the initial state and thus the SSB will not appear. 




FIG. 1: (Color online). The schematic illustration of spontaneous 
symmetry breaking (SSB) in thermalization. (a) The heat bath usu- 
ally does not single out a particular state from the degenerate ground 
states at zero temperature; (b) Breaking the symmetry of the poten- 
tial will cause a preference of the system to a particular ground state 
at zero temperature; (c) After the symmetry is recovered at zero tem- 
perature, the system will remain in the preferred state. 



In this case, the OQS enjoys the so called anti-thermalization 
effect: some information of the initial state is kept after the 
OQS is thermalized to a steady state. 

We will study a three-level system interacting with a heat 
bath. The two lower (or higher) energy states \g\) and \g2) of 
this three-level system are degenerate and can be split by ap- 
plying an external field. The heat bath is modeled as the bath 
of harmonic oscillators. The dynamic process of the system's 
approaching the steady state at zero temperature will be care- 
fully analyzed from the master equation approach. We will 
prove that if there exists a non-vanishing coupling to the bath 
for arbitrary two energy levels of the three-level system, then 
the third law of thermodynamics is valid thanks to the SSB. 
On the other hand, for the conventional A- and V-type atoms, 
we will reveal the exotic anti-thermolization effect. This ef- 
fect happens as a result of the absence of the bath coupling 
induced quantum transition between \g\) and \g2) and the oc- 
currence of the quantum interference between the transition 
to \gi ) and the transition from \g2). At zero temperature, these 
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FIG. 2: (Color online), (a) The schematic diagram of a three-level 
system immersed in heat bath, which consists of a set of harmonic 
oscillators, (b) The energy level diagram of the three-level system, 
arbitrary two energy levels of which couples with the heat bath. 

conclusions coincide with those reached in the context of the 
spontaneous emission of V-type atom in vacuum [8|]. 

SSB in thermalization. — In general, the Hamiltonian Hs 
of an OQS to be thermalized can be written as Hs = 
Y 1 n,a^n\n,a){n,a\, where \n,a) (a = 1,2, ■•• ,d„) are de- 
generate states correspond to the same eigenvalue E„ (n = 
1, 2, 3, ■ • • ) and d„ is the degree of degeneracy. Let an ex- 
ternal field be applied to break the energy level degeneracy 
as E„ — > E„ + A a and then let this degeneracy split sys- 
tem contact with a heat bath of temperature T for a time 
longer than the conventional relaxation time. Then the sys- 
tem is supposed to be thermalized to the thermal equilib- 
rium state p s (A a ,B) = exp[-j8 (E„ + A a )]Z _1 |n, a)(n,a\, 
where B = 1 / {k B T) is the inverse temperature (hereafter, we 
set k B = l,h = 1) and Z = Jr[exp{-BHs)] is the partition 
function of the OQS. 

We observe that for ps (A a ,f3), there exist the following two 
limit processes: 

lim lim p s (A a ,8) = \l,l)(l,l\, (1) 

A„-*0 jS-t+co 

lim lim p s (A a ,B) = - Y \\,a){l,a\, (2) 

/3^+ooA„^0 d\ t—l 

a=l 

where |1, 1) denotes the ground state with vanishing energy 
in the presence of the external field. Note that taking the two 
limits in different orders leads to completely different results, 
the former being a reflection of the SSB phenomenon. 

We would like to remark that though both of the two results 



are correct in the mathematical sense, the former is physically 
more acceptable than the latter, which is in accordance with 
the generic version of the third law of thermodynamics. In- 
deed, in view of the existence of the perturbation breaking 
the degeneracy, it seems reasonable to let the temperature ap- 
proach zero first in the calculation. However, things are not 
so simple. In fact, due to quantum interference effect, even if 
we let the temperature approach zero first in the calculation, 
the happening of SSB is not unconditional. This is the main 
conclusion of this letter. 

Dynamic thermalization of a three-level system. — 
Generally, the process of dynamic thermalization begins from 
a factorized initial state p(0) = ps(0) ® p B (B), where p B (B) 
is the thermal state of the heat bath. Let the heat bath be 
modeled as the harmonic oscillator system with the Hamilto- 
nian H B = YjjWjCL^cij. Then we have p B (B) - exp(-BH B )/Z B 

where Z B = Tr[exp(-BH B )] is the partition function. The 
time evolution of the total system driven by the coupling 
Hi between the system and the heat bath is determined by 
U(t) = exp[-i{Hs + H B + Hj)]. The steady state of the system, 
as t — > oo at T — 0, can then be obtained by calculating the 
reduced density matrix p s (t,J3) = Tr B [U(t)p s (0) ® p B UHt)] 
of the OQS (Tr^ stands for tracing over the heat bath). 

To be specific, let us study the dynamic thermalization of 
a simple three-level system. The Hamiltonian Hs = u)2& e e + 
Aa glgl of the three-level system (as illustrated in Fig. |2|i is 
written in terms of the flip operators & a g = \a)(J3\ (a, ft - 
e >gi,g2), where <±>\ = a> e -a> gl (I = l,2)andA = o>2-cl>\. Here, 
we choose the eigen-energy of the state \g2) as the energy zero 
point. The interaction Hamiltonian of the three-level system 
with its heat bath reads 

Hi-^j °~e gl Bi + d- glg2 B 3 + h.c, (3) 

1=1,2 

where B[ = 2 ; /7/(wj)a 7 (I = 1,2,3). For simplicity, we as- 
sume T]i(a)j) to be real below. 

Under the Born-Markov approximation, the evolution of 
the reduced density matrix of the three-level system is gov- 
erned by the master equation, 
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Ps = -i[p s ,A& g2g2 ] + Yj£iIPs] + £x[psl (4) 

1=1 

where 



-Ci=ia[ps] 
■CslPs] 
■CxlPs] 



■j(n(coi)+ l){2&g ie ps& egl -o- ee ps -psd- ee )+ -^n{ioi){2& eg fi S &g ie - & gl gfi s -p s d- g!gl ), 

-y (n(A)+ l)(2fr g2gl psfr glg2 -fr g]gl ps -p s o- glgl )+ y«(A)(2(7j lS ,p S (J s , Sl - o- g2g2 p s -p s o- g2g2 ), 



^(«( Wl )+l)+?y(»( W 2)+l) 
721 



(d~g,ePSO- eg2 + d- g2e psd- eg] ) + ^y-n(«l) (&e gl Psd- g2 e + &eg 2 PSO- gx e 



- & g i^Ps -Psd- glg2 ) + -^-n{0)2){& egl p s & gie + & eg2 p s & gie - o- glg2 p s -Psd- g2g] ). 



(5) 
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Here, the decay rates yi - 27rg(a>i)\rii(cL>i)\ 2 and y/,„ = 
27rg(a>i)T]i(a>i)ri m (a>i), for W3 = A and Z ^ m, m = 1,2, de- 
pend on the mode density of the heat bath; h(u>) = 
l/[exp(j0w) - 1] is the thermal average excitation number for 
the boson mode of frequency a> at temperature T. Note that in 
the master equation we have neglected the Lamb shifts. 



The evolution of the density matrix elements governed by 
the master equation can be described with the optical 
Bloch equation X = MX, where the state vector X = X(t) 
and the coefficient matrix M are respectively defined as X = 
X s e X s and M = R S, with X s = ((& egl >, {d~e g2 )) T and 



X R = ((n(o>i) + \){cr ee ) - h{(D X ){fr gm >, («(w 2 ) + Y){&ee) ~ n{u> 2 ){& g2g2 ), (n(A) + l)(fr glgl > - n(A)(fr g2g2 ), Re[C], Im[C]) r , 

' -yi(2«(wi) + 1) -y 2 (n(«i) + 1) y^n{o}{) yi 2 n(«i)(n(wi) + 1) + y 2 i«(w 2 )(2n(«i) + 1) 

-yi(n(w 2 ) + 1) -y 2 (2w(« 2 ) + 1) -y 3 n(w 2 ) yi 2 n(wi)(2n(w 2 ) + 1) + y 2 in(w 2 )(«(w 2 ) +1) 

R = n(n(A)+l) -y 2 «(A) -y 3 (2n(A) + 1) yi 2 n(A)n(wi) - y 2 i(n(A) + l)n(w 2 ) 

^ fl 4 4 A 

-A Ass 

-i[ ri (2n( Wl ) + 1) + y 2 (n(co 2 ) + 1) + y 3 (n(A) + 1)] ; 



S = 



-Hit'?) 

ik - ±[y 2 (2n(w 2 ) + 1) + yi(n(wi) + 1) + y 3 «(A)] -^n(^i) 



(6) 



where C = (d~ g2gl ) and /?44 = ^55 = -[y\h(ui) + y 2 n(w 2 ) + 
y3(2«(A) + l)]/2. We can prove M to be negative-definite or 
have vanishing determinant. Thus the optical Bloch equation 
can possess steady state solutions. 

We first consider the thermalization of the three-level sys- 
tem at finite temperature T + 0. In this case, det(M) + 0, 
thus the steady state solution of the optical Bloch equation is 
X = 0, or 

(^g2gi)ss = (d-e gl ) SS = (d-e g2 )ss = (7) 

for I = 1,2. Here (A) ss = Tr s [Ap s (oo)], A being the op- 
erator concerned. Considering the normalization condition 
(d'gigi)ss + (^282)55 + {o~ee)ss = 1, from Eq. (O we then have 

for/= 1,2. 

Next, we consider the case with vanishing bath temperature. 
In this case, we have «(wi) = n(w 2 ) = n(A) = and thus 
det(M) = 0. We thus cannot obtain the steady state solution 
of the optical Bloch equation by simply setting X = as the 
resulted equation = MX would not have a unique solution. 
Instead, we have to turn to obtain the transient solution first 
and then consider the long time behavior. One can expect that 
it cannot be determined without regard to the details of the 
interaction or the initial state. Indeed, complexity will arise 
here. Let us focus on the case with y 3 + in this section, and 
leave the case with y? = to the next section. 

When y 3 + 0, namely, there exists a bath induced cou- 
pling between the states \g\) and |g 2 ), the transient solution, 
which is not presented here for the technicalities, results in 
(d~ee) ss = (°~g lgl )ss = 0, and (d- g2g2 ) ss = 1. This implies that 



whether the three-level system has degenerate ground states or 
not, the steady state of the master equation © is just the ther- 
mal equilibrium state (oc exp(-fiHs)) of the OQS even at zero 
temperature. This conforms to the conventional idea. How- 
ever, when y 3 = 0, in the next section we will see an exotic 
nature of thermalization as T — > 0. 

Now it is easily seen that taking the limits in dif- 
ferent orders leads to the following different results: 
limA^olirrj jg ^ +00 (o- glgI ) M = 0,lim A ^olirrj jg ^ + o ((7- g2 g 2 ) M = 1, 
and Hm^+oolimA-^gigi)™ = hm^+co \im A ^ Q {d- g2g2 ) ss = 
1/2. Following the first procedure, we will reach the con- 
clusion that the final steady state is the pure state |g 2 ) and that 
there exists SSB effect in the thermalization . But if we adopt 
the second procedure, we should conclude that the final steady 
state is the maximally mixing state (\g\)(g\ \ + |g 2 )(,g 2 |)/2. By 
the way we remark that the SSB can also be seen from the 
von Neumann entropy S = — Trs [p5( 00 )lnps( 00 )] of the 
steady state ps(°°) of the three-level system. In Fig. [3] we 
plot the von Neumann entropy S as a function of A and T. 
In the figure the character of double values of S at the point 
(A, T) = (0, 0) is clearly illustrated: along the route (T - 0, 
A — > 0), the von Neumann entropy S — > while along the 
route (A = 0, T — > 0) the von Neumann entropy S — > 1 . 

Anti-thermalization by quantum interference. — We have 
shown that when y 3 + 0, in the dynamic thermalization of 
the three-level system all initial information will finally be 
erased. As mentioned above, things are not so simple when 
y 3 = 0. In fact, the steady state of the A-type three-level 
system immersed in a zero temperature bath will depend on 
its initial state if the bath coupling between the two lower 
levels is forbidden. This phenomenon is referred to as anti- 
thermalization. 

When y3 = 0, the analysis of the A-type three level sys- 
tem immersed in a heat bath can be made in the same way as 
presented above. The time evolution is described by the mas- 
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FIG. 3: (Color online). The von Neumann entropy of the steady state 
density matrix for the three-level system is plotted versus the temper- 
ature T and the energy difference A between the states \g{) and \g 2 ). 
The two red arrows indicate its multi-value feature as both T and A 
approach zero 

ter equation © with 73 = 0. At zero temperature, we have 
det(M) = 0. Thus the optical Bloch equation is reduced to 

(& ee ) = -(71 + j2)(o- e e), (°~g lgl ) = Jlio-ee), 

Re [<^i>] = ^(ri2+r2i)<o- fe }-AIm[<o- S2Sl )], (9) 
Im{& g2gl ) = ARe[<<x g2gl >], (d- egl ) = ~(7i + 72)(d- egl ), 
where I — 1,2. 

These equations can be solved straightforwardly to ob- 
tain the transient solutions for (& ee (f)), (cr glgl (f)) and (& g2g , (f)) 
where I — 1,2. From these transient solutions it follows that 

{*m)» = <<W«)> + «M<>)>, (10) 

7i + 72 

Re[<<W„] = Re[(cr g2g ,(0))] + (y ' 2 + 72l) (& ee (0)\ 

2(7! + 72) 

and (& ee } ss = (& egl ) ss = (d- eg2 ) ss = , Im[(& g2gl ) ss ] = 
Im[(& g2gl (Q))], where / = 1,2. This is just the steady state 
solution of equation (0. 

Equations ( fTTT ) clearly show that the steady state of the 
A-type three-level system depends on its initial state. In 
the steady state, the decaying probabilities from the excited 
state \e) to the ground states \g\) and \g2) are respectively 
7i/(7i + 72) and 72/(71 + 72)- Moreover, as is shown in 
Eq. (fTTl i. in the present case the dynamic thermalization will 
preserve or even increase the off-diagonal elements of the den- 
sity matrix of the initial state while in the previous case with 
73 + 0, the dynamic thermalization will lead the system to a 
steady state whose density matrix possesses no off-diagonal 



elements. It is also noticed that the steady state of the three- 
level system is independent of the initial off-diagonal elements 
between \e) and \gi) (1=1, 2), and when it is initially prepared 
in the superposition state of the two ground states its final 
steady state will be the same as the initial state. Finally let us 
present two simple examples. Take 71 = 72 = 712 =721 =7 
and A = 0, then for the initial state |iA(0)) = \e) we have the 
steady state |iA(°°)) = (\gi) + \g2))l V2~ and for the initial state 
l<A(0)> = \g\) we have the steady state \i//(oo)) = \g\). 

Summary. — In summary, in this letter the SSB effect in dy- 
namic thermalization is studied through a three-level system 
immersed in a heat bath inducing cycle transition couplings. 
Careful calculation is carried out from the master equation 
approach to examine the thermalization dynamics when the 
temperature approaches zero. By this investigation it is con- 
cluded that when there is no selection rule to forbid any one 
of the bath induced cycle transition couplings, the canoni- 
cal thermal state can be reached as a steady state solution of 
the master equation at zero temperature and if the bath in- 
duced transition between the two lower (higher) energy states 
of A-type (V-type) atom is forbidden the anti-thermalization 
phenomenon will happen due to the quantum interference be- 
tween the transition from the lower state and the transition to 
the higher state. In this latter case, the final steady state of the 
three-level system will depend on its initial state, and thus will 
preserve some of the initial information. This means that the 
initial information of the system cannot be completely erased 
by thermalization and the third law of thermodynamics does 
not work in the conventional fashion. 
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